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Abstract 

We discuss the role that higher derivative operators play in field theory orbifold compact- 
ifications on S\jZ% with local and non-local Scherk-Schwarz breaking of supersymmetry. 
Integrating out the bulk fields generates brane-localised higher derivative counterterms to 
the mass of the brane (or zero-mode of the bulk) scalar field, identified with the Higgs field in 
many realistic models. Both Yukawa and gauge interactions are considered and the one-loop 
results found can be used to study the "running" of the scalar field mass with respect to 
the momentum scale in 5D orbifolds. In particular this allows the study of the behaviour 
of the mass under UV scaling of the momentum. The relation between supersymmetry 
breaking and the presence of higher derivative counterterms to the mass of the scalar field 
is investigated. This shows that, regardless of the breaking mechanism, (initial) supersym- 
metry cannot, in general, prevent the emergence of such operators. Some implications for 
phenomenology of the higher derivative operators are also presented. 
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1 Introduction 



The physics of extra dimensions is receiving a strong research interest in the context of effective 
field theory approaches to compactification, sometimes referred to as "field theory orbifolds". 
While string theory may provide the ultimate description of compactification, field theory orb- 
ifolds can also consistently (re)address interesting issues such as supersymmetry breaking, the 
hierarchy problem, radiative corrections or their experimental signatures. In this work we inves- 
tigate the implications of supersymmetry breaking in 5D field theory orbifolds for the one-loop 
corrections to the mass of the 4D scalar fields. 

An aspect that is somewhat overlooked in recent studies of radiative corrections in gauge 
theories on field theory orbifolds (see for example pQ-jHl), is the role of higher derivative oper- 
ators 1 . In general it is thought that such operators, being higher dimensional, are suppressed 
if the scale where they become relevant is high enough. However, from a 4D point of view the 
most natural such scale is the compactification scale 1/R and, if this is low (TeV scale), such 
operators can have a significant effect even at low energies. For generality we shall consider the 
role of such operators in orbifolds with a radius R of an arbitrary value (large or small). 

The models that we investigate are 5D N=l supersymmetric, compactified on S1/Z2 and 
the interactions are localised superpotential (in the extra dimension) and gauge interactions. 
These are standard interactions in any higher dimensional theory which aims to recover after 
compactification and at low energies, the Standard Model (SM) or its supersymmetric versions. 
Such interactions induce loop corrections to the mass of the zero mode scalar field (if this is a 
bulk field) or to that of a brane scalar field. It is important to note that this field is regarded in 
many models as the SM Higgs field EQ-[E]- Therefore such loop corrections are important for the 
hierarchy problem, electroweak symmetry breaking or running of the mass of the scalar field. 

Recent studies of loop corrections to the mass of the 4D scalar field showed that radiative 
corrections are of type (see for example 01 ) 



where rj = —1 (+1) when a is the Yukawa (gauge) coupling. As a result one can have radiative 

electroweak symmetry breaking 2 triggered by towers of Kaluza-Klein states of bulk fields which 

couple, via the aforementioned interactions, to the (zero-mode or brane) scalar field associated 

with the Higgs field. Note that in the absence of a rigorous mechanism to fix 1/R to low values 

(TeV or so) such results do not provide a solution to the hierarchy problem even at one-loop. But 

they open new ways to address these problems and we think this is worth further investigation. 

For works on this topic see for example |10 | -|13 | . 
2 assuming a vanishing mass at the tree level. 
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As we shall see shortly, it turns out that the non-renormalisable character of these theories 
becomes apparent earlier than one may like, but with potentially interesting consequences. That 
means that results like Q are in general altered by the effects of higher derivative operators, in 
some cases already at one loop. Their effect depends on the size of the radius R, but as already 
mentioned, here we keep R arbitrary. The result of our one- loop Yukawa corrections shows that 



where £ is an (unknown) coefficient of the higher derivative operator. 

Eq.(j2J) shows that the scalar mass depends on both R and 1/R which 3 we find rather in- 
teresting. Further, if R is somehow fixed to a large value (inverse TeV scale) to explain a low 
mass for the scalar without large fine tuning, the second term in the first equation becomes 
more important. Given the unknown coefficient £ of the higher derivative operators, this can 
affect the predictive power of the models. Further difficulties at consistency level may arise for 
special values of this parameter, as we discuss in Sectional Conversely, if R is very small (q 2 
fixed), the role of higher derivative operators is suppressed, but then the first term re-introduces 
the quadratic mass scale (hierarchy) problem at one-loop, familiar from the Standard Model. 
While this is the general picture, a detail analysis should consider the 0(q 2 R 2 ) terms in iffi) 
whose expression can be computed in general cases using our technical results. Finally, our 
calculation can be used to re-address previous studies of the radiative electroweak symmetry 
breaking induced by towers of Kaluza-Klein modes, and to investigate the one- loop running of 
the scalar mass and its ultraviolet behaviour under the UV scaling of the momenta, q 2 — > aq 2 . 

We investigate whether the presence of higher derivative operators and the above consider- 
ations depend on the way 5D N=l supersymmetry is broken. This is the main purpose of the 
paper. We consider both local and non-local supersymmetry breaking, transmitted to the visible 
sector via radiative corrections. The non-local breaking includes discrete and continuous Scherk- 
Schwarz twisted boundary conditions for the bulk fields. The plan of the paper is as follows: 
in Section El we discuss the orbifold compactifications and supersymmetry breaking by local F 
terms of a brane field at a hidden brane (Section I2.1j) and non-local breaking (Section 12.2(1 by 
Scherk-Schwarz boundary conditions. We show that higher derivative operators are generated 
as counterterms for Yukawa but not for gauge interactions. It turns out that the presence of 
these counterterms is rather independent of the above ways of supersymmetry breaking and this 
is discussed in Section |2 Some phenomenological implications of the higher derivative operators 
that are found are discussed in Section |IJ Appendix contains formulae relevant for the study 
of the running of the scalar field mass (with the momentum scale), which may also be used in 
other applications. Appendix El outlines a dimensional analysis needed in Section^ 

3 Such dependence of loop corrections on both R and 1/R is familiar in one-loop string corrections (T duality) 
such as those to gauge couplings |14j . 
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2 Higher derivative counterterms on S1/Z2. 



On the S1/Z2 orbifold, one can consider vector supermultiplets and hypermultiplets. The former 
may be described in a 4D language as made of vector superfield V{\\, A^) and an adjoint chiral 
superfield S((<r + iA§)/-\/2, A2), where Ai^ are Weyl fermions, a a real scalar and A^,A$ the 
5D gauge field. The hypermultiplet contains two chiral superfields VO an d & c ((j) c ,ip c ) with 
opposite quantum numbers, and where 4>, 4> c are complex scalars and ip, ip c are Weyl fermions. 
The orbifold conditions considered are such as the gauge field A^ has even parity (has a massless 
zero mode) to respect the 4D gauge invariance. We consider the following parity assignments 

®(x,-y) = ®(x,y), V(x, -y) = V(x, y) 

$ c (x, -y) = -$ c (x, y) £(x, -y) = -E(x, y), (3) 

where $ can be any of the SM fields Q, U, D, L, E. As a result, the original 5D N=l super- 
symmetry is broken and the fixed points {y = 0, 11R) of the orbifold have a remaining 4D N=l 
super symmetry. Further, we consider the following localised interaction, 



£4 = JdyS(y)l-Jd 2 e [\ t QUH u + X b QDH d + ---] +h.c.|. (4) 

The 5D coupling X t = f B ,t/M? = (2vri?) n / 4ii where f 5>t (/ 4 , t ) is the dimensionless 5D (4D), M* 
is the cutoff of the theory. In the following Q, U, D superfields are always bulk fields and have 
mass dimension [Q] = [U] = [D] = 3/2, while the Higgs field H u d can be a brane field [H u ^ = 1 if 
n = 1 or a bulk field 4 = 3/2, when n = 3/2 (when it also has a H c a d partner). The above 

spectrum and interactions define our minimal model. While this is not phenomenologically 
viable, it is sufficient to illustrate the idea of the paper. Moreover the spectrum and interaction 
considered are generic in many detailed models which reproduce the SM or its supersymmetric 
versions. Such models, for which our findings are relevant, can be found in refs.pQ to [H]. 

The purpose of this work is to investigate the correction to the mass of the brane (or zero- 
mode of the bulk) scalar fields 4>h u , <fe d induced 5 by towers of Kaluza-Klein states of Q,U,D 
superfields, via interaction @. Gauge corrections will also be considered. Since similar consid- 
erations apply to both H u and Hj, we shall compute the one- loop correction to the mass m ( f >H of 
4>h u , hereafter denoted simply 4>h- To compute the one-loop corrected m$ H one must first eval- 
uate the spectrum of Kaluza-Klein modes, using eq.©. However this spectrum also depends on 
the mechanism of the further breaking of 4D N=l supersymmetry which is then transmitted to 
the "visible" sector at y = 0. We shall distinguish two cases discussed separately: I. Localised 
supersymmetry breaking and II. Non-local supersymmetry breaking. 
4 If H Ui d are bulk fields they also satisfy a condition similar to that for $ in eq.Q- 

5 One needs two Higgs fields of opposite hypercharge to avoid (quadratically divergent) FI corrections |15) . 



4 



2.1 Higher derivative operators from localised supersymmetry breaking. 

We consider that supersymmetry is broken at a distant (hidden) brane located at y = irR by 



£ 4 = J dy 5{y - ttR) J d 2 9 Ml Z + h.c. 



(5) 



where Z is a (gauge singlet) brane field at y = ttR. The bulk fields Q, U feel the supersymmetry 
breaking via the couplings 



dy5(y - ttR) 



dH 



ZLQlQZtZ + ^LrtUZtZ 



(6) 



Therefore, when (Z) ~ Fz6 2 , the bulk fields such as (f>Q } u with non-zero coupling at the y = ttR 
brane have the spectrum modified. One can show that (for details see for example 3 ) 



tanf-7r Rm 



4>M,n\ 



CM Fl 



in 



4>M,n 



( n+ 2)R 



M = Q,U. 



(7) 



2 M> fa ; 

With the choice cm ~ Fz ~ Mi? > an d to leading order in l/(RM*), one can set u = 1. This 

is usually referred to as "strong" supersymmetry breaking, otherwise u is a series in l/(i?M !)! ). 

Unlike the fields 4*q,u-i their fermionic partners ipQ t u do not couple to the vev of Z, and their 
mass is not changed by ©. This results in the breaking at y = ttR of the remaining 4D N=l 
supersymmetry of zero modes. Similarly, the fields tpg v , 4>q v do not couple to Z (see eq.©), 
thus their mass spectrum is not affected. Finally, the Higgs field H u (also H^) is considered to 
be a brane field, localised at y = 0, so it does not have tree level interactions with the field Z. 
However, it feels the supersymmetry breaking at y = irR through loops of the bulk fields Q, U, 
(or Q,D for H4) via eq.®. Therefore the spectrum of the bulk fields is, using Q, 



k_ 



u 



R 



k > 0, 



k > 0, 



R> "Z 1 



M = Q,U. (8) 



where = 1/2 and u = 1 for strong supersymmetry breaking. The case with u=f\ is discussed 
in Section |31 Finally, the wavefunction normalisation coefficients (at y = 0) are rft M = 1, rft u = 
r]^ M = l/\/2 (A; > 0). These results will be used in Section \2. 1.1 1 to compute fn<f> H at one-loop. 

Finally, let us also consider the case of supersymmetry breaking at y = ttR with a brane- 
localised gaugino mass on S 1 /Z2, 

1 



C 4 



dy 5{y - ttR) / d z 6 



c w 



ZTr 



W a W n 



+ h.c. 



(9) 



(4ff 5 ) 2 Ml 

Here W a contains gaugino as its lowest component. Note that only the even parity gaugino Ai 
(see eq.©) couples to Z and feels the supersymmetry breaking at y = irR. Further, Ai and A2 
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are coupled through their kinetic term. Eq.@ is the counterpart of eq.© for the gauge case. 
Using the equation of motion one finds [H] 

c w 



tan[7ri?m> 



4M2 



Fz, 



Then, the spectrum of gaugino is 



with 



k + P 1-^7 

m Xl = m\ 2 = — — , k G Z, 



1 . (cwFz\ 
p = — arctan . 



(10) 



(11) 



(12) 



This result is used in Section [2.1.21 for the one-loop gauge correction to the scalar mass m ( j >H . 

2.1.1 One-loop mass correction due to Yukawa interaction 

In the on-shell formulation, the interaction in eq.@ becomes in component fields [3] 



EE ^h, 

k=l 1=0 

oo oo oo 



m 4> c Q ,k rik Qr )f U <pQ,k </>u,l 4>H + h.c. + (Q<r->U) 



- E E E ( 2 A*) 2 hN*(^)VU*w*Wfl + (q - v) 



k=0 1=0 m=0 

oo oo 



- E E ( 2 At) [Vk^^kMH + h.c. 

fc=0Z=0 



(13) 



Note that the sum over k in the first line is from k > 1 since the field </> c is odd under orbifolding, 
eq.(j3J). If we work on-shell (off-shell) there are three (two) types of diagrams contributing to the 
mass of brane Higgs field (pg , see Fig. ^ The one- loop contributions in Euclidean space are 6 



B 



i{2U,t?N c Y, [ 



Q„4>v 



k>0,l>0 



d d p 



1)0 + qf I* 



2 ,,4-d 



(27r) d ((p + g) 2 + m2 )(p2 + m 2^, 



(2/ 4 ,t) 2 iV c E [ 



fc>0,Z>0 



2p.(p + <?) //■ 



4-eZ 



(2^ ((p + q) 2 + ml Q>k )(p 2 + m 2 M ) 



+ (Q^U) 



(14) 



This form of the radiative corrections can be shown to be equivalent to that derived using 5D Green functions 
in mixed positions-momentum space (computed in ref. ,3 , eqs.(57), (60), (64)), and evaluated at y = 0. This also 
suggests the (brane) localisation of the corresponding counterterm, see later. 
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F ¥ 

r Q,k r U,k 

Figure 1: Diagrams contributing to the two-point Green function of the (brane or zero mode of the bulk) 
scalar field <f>jj- For the left diagram one should also add the similar contribution with Q «-> U. 

Here d = 4— e (e— »0), N c is the number of colours, fi is the scale introduced by the DR scheme, 
q 2 is the external momentum and the two double sums extend to infinity. The index B, (F) 
stands for bosonic (fermionic) contributions. One then uses the spectrum (JHJ) and coefficients r\ 
given after this equation, perform the integrals over p in DR, then the double sums, to find 



(2/ 4>t ) 2 iV c 



1 f d d p 



fc>0,Z>0 



(-2)(p + q) 2 i& 



2 S k ,oJ (2ir) d ((p + qf + k 2 /R 2 )(p 2 + (I + l/2) 2 /R 2 ) 



(2/ 4 
2(4 



(2U, t ) 2 N c 



^^N c £dx{ J 2 [l/2,0,c]+q 2 R 2 (l-x) 2 J 1 [1/2, 0, c] } 

1 1 f d d p 2p.(p + q) fi A ~ d 



fc>0,Z>0 



2 s w 2^.o J (2vr) d ((p + q) 2 + k 2 /R 2 )(p 2 + / 2 /i? 2 ) 



(2/ 4 ,t) 2 «. 



2 (47T.R) 



2 iV c / da; 



2-e/2 



7T 



J 2 [0,0,c] +g 2 i? 2 x(x- 1) Ji[0,0,c] }, (15) 



where K e = (2TTfiR) € . The functions J12 have the following definition and leading behaviour in e 



Jj[ci,c 2 ,c] = J2 
ai = 



(ft 



v-*l 2 

1-x), 



- 7r t (c+ai ( fei +ci ) 2 +a 2 ( fc 2 +C2 ) 2 ) 



(-vrc) j 



J roi 



a 2 = x, c = x(l — x) q 2 R 2 . 



+ O(e ), j = l,2. 



(16) 



For a complete expression of the functions see the Appendix, eas. ()A-l|) to (|A-3|) and (|A-11|) 
to (|A-15jl . It is important to notice that the leading (divergent) behaviour of J\$, depends on c 
and a\,a 2 but is independent of ci, c 2 - The dependence on c is very important, since it is only 
for c ~ q 2 R 2 7^ i.e. non-zero external momentum that one is able to "see" the poles of Jt,2- 
After adding the bosonic and fermionic contributions in l)15[). one finds 
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m: 



2(4^) 2 c 



dx(2/n) J 2 [0,0,c] -J 2 [1/2,0, c] 



x(x-l)Ji[0,0,c] -(1 



: Ji [1/2,0, c] 



(17) 



with c as in (|16|). Note that if q 2 = the second line above is absent, so m^^q 2 = 0) is given 
by the first line alone. Further, in the difference ^[0,0, c] — ^[1/2,0, c] the divergent part 
q 4 R 2 /e in each J2 cancels away to give a one-loop finite m^^O) ~ 1/R 2 . This cancellation is 
ensured by the equal number of bosonic and fermionic degrees of freedom, enforced by the initial 
supersymmetry. Using the leading behaviour of the functions J7i,2 we find 



At 



(2/ 4 ,t) ; 
2 8 



N c (q 4 R 2 ) 



\n(2nRn) 



+ j~ 2 Oiq'R 2 ) 



(18) 



where the /^-dependent term in the square bracket shows the regularisation scheme dependence 
induced by K e . The finite part 0{q 2 R 2 ) can be evaluated from the second line in (|T7|) using 
eas. (|A-2|) . (|A-3|) in the Appendix to give the full running of the mass wrt momentum scale q 2 . 

Eq. (jl8j) shows the presence in the sum of bosonic and fermionic contributions, of a pole 
multiplied by quartic dependence on (external momentum) q, originating from the two J7i's. 
The result is that the one loop scalar mass is not finite and has a UV divergence similar to that 
cancelled by (initial) supersymmetry in the J 2 dependent part. Therefore one must add in the 
action a higher derivative counterterm 7 to m 2 ^ H (recalling that (j)H is a brane field) 



(19) 



<Fxd z 6d A 9\i W U UH U ~/£ t d*xR z ct) l H n z (j) H + ... 



The presence of the UV divergence and of corresponding higher derivative counterterm shows 
that, although the initial theory was supersymmetric, its non-renormalisable character is never- 
theless manifest through a counterterm generated by the large number (multiplicity) of Kaluza- 
Klein modes which contribute to m ( / )H . 

In Table^we provided a general dimensional analysis for when higher derivative counterterms 
can arise from a localised superpotential. The table shows the number of loops in perturbation 
theory at which such counterterms can be generated, by assuming that any of the fields in the 
superpotential can be bulk or brane fields. This information can be used in realistic orbifold 
models, to avoid such operators at low orders in perturbation theory. 

7 This result is similar to that in llol which had instead 5D N=l supersymmetry broken to N=0 on S\ j {ZiV.Z 1 -}) . 
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[Q] 


[U] 




[A] 


# (no. of loops) 


3/2 (bulk) 


3/2 (bulk) 


1 (brane) 


-1 


n = 1 


3/2 


1 


3/2 


-1 


n = 1 


3/2 


3/2 


3/2 


-3/2 


n = 1 


3/2 


1 


1 


-1/2 


n = 2 


1 


1 


3/2 


-1/2 


n = 3 



Table 1: This is an estimate of the number of loops n when higher derivative operators may be generated. 
The localised superpotential J dy §{y) d 2 9 XQUH U with Q,U,H U as brane/bulk fields generates higher 
derivative counterterms to the scalar zero mode of H u (if a bulk field) or scalar component of H u (if a 
brane field). An example is J d 4 d {\ 2 ) n H u UH u . The table is a dimensional estimate of the number of 
loops when this counterterm arises located at y = 0, in function of the nature (bulk/brane) of the fields. 



2.1.2 One-loop gauge correction to a brane scalar mass 



Let us now consider the one-loop gauge contribution to the mass ?7U„ of the brane scalar §n 
located at y = 0, assumed to be charged under a U(l) group. This is induced by the action (jHJ) 
and in the following we use eas. ()lU|) to ()12|) . In the dimensional regularisation with d = 4 — e, 
bosonic (gauge) and fermionic (gauginos) contributions to the scalar self-energy at nonzero 
external momentum q 2 , are respectively 



n.ez J 



d d p 



p.(q+p) 



(27r) d (p 2 +n 2 /R 2 )(q+p) 2 



and 



E 



d d p 



p.(q + p) 



n __ z ., (27ry(p 2 + (n + p) 2 /R 2 )(q + py 
Then, we find the one- loop corrections as 



(20) 



(21) 



2 / 2n 

m <p H {q 



and 



b = 9 ^§- r dx [( 2 - \) ^ ° ] - ™ (i - x) q2R2 ^ c] 



(22) 



(23) 



with c = x(l — x)q R and we introduced 
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Gj\pA = 2 / 



dt 



e -^[c+/?(n+p) 2 ] ) j = l )2 ; /? = X. 



(24) 



After some calculations one obtains [HU] 



G2[P,C] 



In 



2sin7r(p + i\J c/ (3) 



e 2m(p+Vc//3)^ +c c I, ( 25 ) 



where c, /? > and Li CT (x) = J2k>i xk ■ Although these expressions are finite (no poles in 
e), whenever one removes a mode from the series, poles arise and this justifies keeping the e 
dependence explicit in their definition ea.(|24jl. 

Therefore, the resulting one-loop correction for the brane scalar is given by 



m l H (<7 2 ) = m% B (q 2 ) b + m\ H (q 2 ) 



f 2tt 3 R 2 J 



dx 



G2[0,c] - G2[p,c] 



4tt 2 R 2 



(q 2 R 2 ) I dxx(l-x) 
Jo 



Qi [0, c] - Qi [p, c] 



(26) 



This result can be further simplified, but for our purpose it is enough to notice that it is finite 
(has no poles in e). Therefore we conclude that no higher derivative counterterms are generated 
at the one-loop level. We note that the momentum-independent mass correction is given by 



4vr 4 i? 2 



C[3]-i(Li 3 (e 2 ^)+ C .c. 



(27) 



This result agrees with the one obtained by changing the infinite Kaluza-Klein sum into a contour 
integral [T7j . 

Alternatively, one can investigate the gauge correction by considering supergraphs. Given 
the presence of subtle differences between component and supergraph formalisms 8 we now use 
the supergraph formalism to show again that no higher derivative counterterms are present. 
This will confirm our conclusion obtained in the component formalism. 

For this purpose we compute the gauge correction to the propagator of a (massless) brane 
chiral multiplet H in the absence of supersymmetry breaking. To do so we need to consider 
only one supergraph with brane-chiral and bulk- vector multiplets "running" in the loop |18j . 
We assume that, as in the 4D case, the soft breaking does not renormalise the propagator of a 



related to gauge fixing in the WZ gauge in component formalism as compared to superfield gauge fixing in 5D. 
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massless brane chiral multiplet. With an appropriate gauge fixing term |19j . i.e. the 5d version 
of the super Feynman gauge, the action for the vector superfield is 



(28) 



£ 5 = / d b xd z 9d 2 9V[-a - d£}V. 
Thus, the propagator for the bulk vector multiplet on S 1 jZi satisfies 

(-q 2 + d!)A 5 (y -y',6- 9') = — ~ ^ ~ 2lTnR ) ~ 

1 n&Z 

Therefore, for < y — y' < ttR, the mixed position-momentum propagator is given by 



(29) 



A 5 = --G 5 (y,y')5\e-e') 



with 



G 5 (y,y') 



cosh[g (y — y' — 7ri?)] 



2 q sinh(-7ri? q) 

In particular, the bosonic part of the propagator at the origin with y = y' = is given by 



G 5 (0,0) 



1 



-E 



2q tanh(vr.R q) 2ttR q 2 + (n/R) 2 ' 



(30) 



(31) 



We obtain the one-loop gauge correction to the propagator of the brane chiral multiplet located 
at y = as 



d^q 



where 



A(q) 



2ttRJ (2tt) 4 
d d k 



A(q) 



n€Z 



d A 9H(-q,9)H{q,9) 



1 



(32) 



(2ir) d {q + k) 2 {k 2 +n 2 /R 2 ) 
v dt 



(ATT) d / 2 (7Tfl 2 R 2 ) d / 2 - 



t d/2-l 



-irt(c+xn 2 ) 



(33) 



where c = x(l — x) q 2 R 2 . Finally, with d = 4 — e and using eqs. (|24j) . (|25|l. one finds that the 
one-loop gauge correction to the kinetic term of the brane scalar (f>H (component of H) in 1)32(1 
has a momentum dependence of type 



q 2 A(q) 



1 1 

(4tt) 2 (irRf 



4 |T 



C[3] + 2yLi 2 (e-^)+Li 3 (e-^) 



+ C(e), y = irRjq 2 . (34) 



Since the above result has no poles in e, we find, using this time the supergraph computation, 
that the wave function of the brane chiral multiplet is not renormalised. Therefore no higher 
derivative counterterms arise at one-loop order, in agreement with the previous computation 
using the component field formalism. 
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Finally, let us discuss the possible higher dimensional (derivative) operators which can be 
induced by the gauge corrections at higher loops, along the lines discussed in Table ^ The 
corresponding operator for a brane chiral multiplet is 



J d A 6{g 2 ) n HUH 



with n being the number of loops. Then, given the (mass) dimension of the higher dimensional 
gauge coupling, i.e. = — 1 in 5D higher dimensional operators can in principle be generated 
for n = 2 (two-loop). However, for a 6D case, [g§] = —2, higher derivative operators may in 
principle be generated for n = 1 (one-loop) 9 . 

2.2 Higher derivative operators from non-local Super symmetry breaking. 

So far we have considered a brane-localised supersymmetry breaking, eqs.©, ©> ©■ I n the 
following we consider a non-local supersymmetry breaking mechanism, with: (1) discrete and 
(2) continuous twisted boundary conditions, which we examine separately. 

2.2.1 Discrete Scherk-Schwarz twists. 

First, let us impose that the 5D fields acquire under a 2ttR shift, a phase which is the i?-parity 
charge of these fields. The action of the i?-parity operator is 

Z 2 ,rQ(x, y, 9) = -Q{x, y, -9), Z 2 , R Q c (x, y, 9) = -Q c (x, y, -9), 

Z 2 , R H(x, y, 9) = H(x, y, -9), Z 2 , R H c (x, y, 9) = H c {x, y, -9), 

Z 2jR V(x,y,9) = V(x,y,-e), Z 2jR X(x, y,9) = E(x,y, -6) (35) 

One also has a condition for U, U c superfields similar to that for Q, Q°. The condition for H, H c 
stands for both H u ^H^ d and applies only if these fields are bulk fields. Eas, ()35|) . © give the 
spectrum relevant for our purpose 

k k 
m 4>M,k = k > 0, m r M ,k = jj, k>l 

m Huk = fc+ R 1/2 , k>0, = k+ ^ 2 , k>0; M = Q,U. (36) 

and rj^ M = rjf M = 1 and r]f M = l/v / 2 <5fc '°. Also if the Higgs field is a bulk field, m^ Hk = k/R 

{k > 0) and m 0HC k = k/R (k > 1). Finally m A ^ k = k/R, m Xl 2 h = (k + 1/2) /R (k > 0). 

Note that the spectrum on the orbifold S\/Z 2 with the R-parity ()35|) has similarities with 

that in the case of S\/{Z 2 x Z' 2 ) orbifold (see for example [3]) with the Z 2 identified with a 

9 A component field computation on T 2 /Z2 shows that one-loop gauge correction to the self-energy of a brane 
scalar is finite |20) . The agreement of this result with that using a supergraph approach is studied elsewhere 
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Z2 i r Scherk-Schwarz breaking of super symmetry. As a result the one-loop corrected m ( p H (q 2 ) 
is expected to be similar to that in S\/(Z2 X Z 2 ) studied in 10 13 j. This similarity is only 
present when considering interactions from one fixed point, and breaks down when overlapping 
interactions from different fixed points are included, as can happen at higher loops. 

The action is similar to that in eq. (|13j) with the remark that the first sum over k starts from 
k = 0. With this information we can compute the one-loop corrections to the mass of the scalar 
component (f)jj of H u , or of its zero mode if H u is a bulk field. One has in this case a result 
similar to eq.(|14|). but note that the wavefunction coefficients tj Fm and the spectrum (m^) have 
changed. In Euclidean space 



im l H (q 2 ) 



B 



;/t 2 ^c£ [r 



fc>0,/>0 
k>0,l>0 



2 f d d p 



2 (p + q) 2 H 



2 ,A~d 



(2^Y{(p + q ) 2 + ml hk )(p 2 + ml U:l ) 

2 r d d p 2p.(p + q) 

(2tt) 4 ((p + q) 2 + m\ Q k \p 2 + m^) 



, (37) 



with the notation ft = 2 n /4t, where n = 1 (n = 3/2) is H u is a brane (bulk) field 11 . In the 
first equation we used that the masses of (j)Q, 4>q, <pu-> 4>u are au equal, and this explains the 
presence of a factor 2 in the numerator of the integrand. 

The calculation of ea, ()37|) with replacements (|3*o]l . proceeds as in Section l2.1.1l fsee also [T3*]). 
After adding the bosonic and fermionic contributions, one obtains 



f 2 

H N c / dx 



16vr 3 ii 2 



o 



J 2 [0,0,c] - J 2 [l/2,l/2,c] 



ft 



+ ^r^? K £ N c<i • dx 

32 7T Z 



x(x- l)Ji [0,0, c] - (1 -x) 2 Ji[l/2,l/2,c] (38) 



The terms involving Jj [1/2, 1/2, c] , j = 1,2 account for the bosonic contribution, while Jj [0, 0, c] 
account for the fermionic part. The functions are given in eq. (|16j) . Comparing (|38|) to (|17j) 
one notices a similar structure, but there is a difference in the arguments of J7i2 hi the two 
equations. Since the pole structure of J\^\c\i C2, c] does not depend on the arguments c±, C2, the 
discussion of the UV divergences is not changed from that of eq. (|17j) . Eqs.(|38|). (|T6|) give again 



A D 2 



+ ]n(2irRfj.) 



+ ^0(q 2 R 2 ) 



(39) 



For the completeness of our analysis we include this case here in detail. 
l If H u is a bulk field, the above correction refers to its zero mode scalar. 
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where 0(q 2 R 2 ) terms are due to J\ functions and can be evaluated numerically using the full 
expression of J\ given in the Appendix. The above result looks similar to that in (|18[) . but the 
exact expression of rrii (q 2 ) is different due to J7i's of different arguments. The counterterm 
has then the structure (when H u is a bulk field) 



Jd A xdyJd 2 ed 2 9 5(y)\ 2 HluH u ~ fffd^xR 2 ^ 4,n a V P 

n,p>0 

~ f 2 Jd 4 xR 2 ^ Hfi n 2 Hfi + ... (40) 
with [At] = —3/2. If H u is a brane field instead ([H u ] = 1 and [At] = —1), the counterterm reads 
J d 4 xd 2 6d 2 8 X 2 HluH u ~f?Jd*xR 2 <l>\ I n 2 <l> H + ... (41) 

2.2.2 Continuous Scherk-Schwarz twists. 

Instead of a discrete twist eq. (|35[) , in this Section we impose continuous twists on bulk fields by 
using the SU(2)r global symmetry. The SU (2)r action under y — > y + 2irR is 

( X ^(x,y + 2nR) = ^ ( M (^), (42) 
A N (x,y + 2irR) = A N (x,y), N = fi,5. (43) 

(+»yx, y + 2«R) = e-^^y X ,y), (44) 

(^)(x,y + 2irR) = (*%)(z tV ), M^Q,U. (45) 

where ((f) m, ^m) anci (^M, V'm) are bulk quark multiplets, M = Q, U. If we also allow the Higgs 
multiplet(s) to live in the bulk, its boundary condition is 

[*%j(x,y + 2nR) = ( | W (46) 

^5 j( X ,y + 27rR) = -(pA(x,y) (47) 

where we note that the higgsinos acquire only a phase of -R-parity because it is a singlet under 
SU (2)r. (If there are two Higgs multiplets in the bulk, one can also use SU(2)h flavor symmetry 
to impose boundary conditions 5 ). 
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The squarks (also the Higgs scalars if they are bulk fields) with a continuous Scherk-Schwarz 
phase have the mode expansion given by 

{%)^y) = ^tjn{y)Vn{ X ) (48) 

where 

Un = e -iu, ny /R f co S (ny/R) \ 
y s'm(ny/R) J 

and (□ — M 2 )(p n (x) = with M 2 = (n + u>) 2 / R 2 . Here we note that the orthogonality is defined 
for eigenstates of SU(2)r doublet as 

1 r^R 

dy {u n {y)) ] u m (y) = S nm . (50) 



2nR jo 

The spectrum of the bulk fields is 

k k 
m V>M,fc = ft> k>0, m 1 pc i k = — } k > 1, (51) 

m<t> M ,k = m ^ M ,k = ^-, k£Z, (M = Q,U). (52) 

k + 1/2 

mil>B,k = = — ^ — , k > 0, (53) 

m<f> H ,k = rn^ H)k = ~ , k e Z, (54) 

and rj^ = rjf = l/y/2 and rjf = \j\ft kfi . Finally, for gauginos one has that m\ 1 2 ^ = (k + tu)/R, 
(k € Z). We note that the zero mode of a bulk Higgs scalar acquires a tree-level mass of u/R. In 
this case, the one-loop Yukawa correction must be larger than this tree-level mass for electroweak 
symmetry breaking 0. However, if the Higgs multiplet is a brane field, such situation is avoided, 
and one can still assume that there is no tree-level Higgs mass. 

The action is in this case similar to that in ea. (|13|) except that the sums over k and I should 
be taken over the whole set of integer numbers. Then the one-loop mass corrections are 

\2 ,,4-d 



-im\ H (q 2 ) 



b 1 1 \hz [Vk 111 1 J ( 2 -) 4 «p + 1) 2 + m i h , k )(p 2 + ' 

-if 2 N V [nVl 2 / ^ 2p.(p + q)^- d 
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Adding the bosonic and fermionic contributions, one obtains the one-loop mass correction as 



f 2 

il AT 

16vr 3 i? 2 c 



dx 



J 2 [0,0,c] -J 2 [u,u,c] 



+ 



f! 



32 n 2 



K e N r 



dx 



x(z-l)Ji[0,0,c] - (1 -x) 2 Ji[(j,a;,c] 



Therefore we find again 

mJ H ( g 3 )=mJ H (0) 



§N c q*R 2 



1 



+ ln(27ri?^) 



(56) 



(57) 



where 0(q 2 R 2 ) terms are e-independent, are due to the two J\ functions in (|o^|) and can be 
evaluated numerically using eas. ()A-2|) . (|A-3|) . Higher derivative counterterms are again required, 
and the same arguments as in eq. (jlU|) , (jlT|) apply. 

In this section we considered so far Yukawa corrections only. Regarding the gauge correction, 
in both cases of discrete and continuous Scherk-Schwarz twists, the spectrum and brane coupling 
of gaugino are the same as in the local supersymmetry breaking. Hence the resulting one-loop 
correction to a brane scalar mass is the same ^7j as before, eqs. ()26|) and (|27|) . 



3 Further remarks on higher derivative counterterms. 

In this section we discuss further the origin of the higher derivative counterterms and their 
relation to the local and non-local character of supersymmetry breaking of Sections 12.11 12.21 

Let us recall first that the radiative correction to the scalar mass from the Kaluza-Klein 
modes has the general structure given by (|58|) below. This equation is just a generalisation of 
eqs. (jTi|) (|T5^ . (|37j) . for m ( f )H , all recovered for particular values of wavefunction coefficients 
a\ : 2 = 1,2 and of "mass shifts" c\ 2 - After a long calculation one obtains 12 (d = 4 — e) 



= E 

fci>0;fc 2 >0 

1 2 



(4tt) 2 e 









<5fc 2 ,o i d<i p 




ap 2 + p(p.q) + 7<? 2 + $ 


_ 2 




_ 2 _ 


J (27T) d 


[{q+p) 2 


+ (k 2 + c 2 ) 2 /R 2 ] [p 2 + (ki + Cl ) 2 u 2 /R 2 ) 


{( 


5 + q 2 ^ 


- 0/2)) ( Cl - 


o\ — \ \ 
2 J 





a u 



6R 2 



cr 2 — 1 



UCl \ C 2 

^-5q 2 R 2 + ^-(a-AP + 8l )q'R 2 
62u 2 a u 



(l + 3ci(l - CTl) + 2cfj + (ci ^C 2 , <71 ^£J 2 , u- 



+ O(e ). 



2 To derive ea. l|58^ one can use the method outlined in Appendix C of ref.|13|. 



b 

u ' 



(58) 
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This result is computed in the DR scheme and as usual, e is a regulator of both the integral and 
the double series in front of it; e is thus a genuine 6D regulator rather than a 4D one 13 . One 
firstly performs the momentum integral to obtain a double series whose summand (function of 
k\,k2) has powers involving e and is multiplied by Gamma functions of e-dependent argument. 
The series are analytically continued and their leading contribution to 0(e ) is obtained. Finally 
one takes account of the Gamma functions and finds the above result 14 . Note that the sums in 
are restricted to positive integers only, unlike the expressions of used previously and 
which involve double sums over Z. The motivation and the advantage of using eq.(|58|l is that 
(unlike the analysis using J1.2) it will allow us to see explicitly how quadratic divergences cancel. 

The divergences in (|58|) are: q 2 /e which account for wave function renormalisation, the terms 
l/(R 2 e) which account for quadratic divergences and finally, the terms q A R 2 /e which account 
for quartic divergences (higher derivative counterterms). 

We can now apply the result (f55|) to the calculations in Sections 2.1 and 2.2 to discuss 
the origin of the higher derivative operators and of other divergences present. For the case in 
Section [2.1.11 for the contribution in the first line of eas. 1)14(1 . (|15|) one has 

Bosonic part : o\ = 2, 02 = 1, ft = 2, 7 = 1, c\ = 1/2, C2 = 0. 

Fermionic part : a\ = 1, cr 2 = 1, f3 = 1, 7 = 0, c\ = 0, C2 = 0. (59) 

while for the case in Section 12.2.11 eq. (|H7j) one has 

Bosonic part : u\ = 2, 0% = 2, /3 = 2, 7 = 1, c\ = 1/2, C2 = 1/2. 

Fermionic part : u\ = 1, 02 = 1, /? = 1, 7 = 0, c\ = 0, C2 = 0. (60) 

In both cases a = 1, 11 = 1, 8 = 0. 

Finally, for the case in Section 12.2.21 with a continuous Scherk-Schwarz phase, one has the 
same expression for the fermionic part as in (|6U|) . but the bosonic part is given by 



£[q 2 ; ci = c 2 = tSj + £ (g 2 ; c\ = u, c 2 = 1 
+ £(V;ci = 1 — o»,c 2 =o>) +£(q 2 ;ci 



u> 



C 2 = 1 — U 



(61) 



where o"i = o"2 = 2,/3 = 2,7 = l and a = l,u=l,5 = in each term. Each contribution £ in (J61J) has 
quadratic divergences l/(R 2 e), as seen from eq. ([58|) . However they cancel in the sum of the first 
two terms (of arguments C2 = oj and C2 = 1 — uj) and also in the sum of the last two terms. This 
is essentially due to summing over the whole set Z of modes in ea.(|48|). Given the above values 

13 If the series in ea. l|58|l were restricted to a finite number of modes then e would act as a 4D regulator only. 
14 In a similar way one can also find the finite, O(e ) terms of l|58|l . 
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of 01,2,01,2, for all cases considered in ea. (|59j) . (JHOJ), (foTj) . one concludes that the coefficient of 
the quadratic divergences l/(-R 2 e) vanishes separately for the bosons and fermions. Note that 
the values of Cj = 0, 1/2 are special for in that case (with corresponding values of <7j), £ alone 
has no quadratic divergence. 

For the origin of higher derivative operators it is important to notice that the coefficient of 
q 4 R 2 /e term is independent of ci,2 and also of o~\ 2. Note that ci,2 which enter the mass formulae 
for the Kaluza-Klein states are in fact set by the boundary conditions for the hypermultiplets 
with respect to the compact dimension. Therefore, this coefficient is independent of the phase 
that hypermultiplet fields have with respect to this dimension and, to some extent 15 , on the way 
supersymmetry is broken. The coefficient of q 4 R 2 /e depends only on a, f3, 7, which in turn 
are controlled by the nature (fermionic/bosonic) of the component fields (via their propagator 
in momentum space). Technically, the term q 4 R 2 /e is strictly the result of the presence of two 
sums in front of the integral (|58|) . over terms with ki 7^ 0, i = 1,2. (this explains the absence 
of such divergences and corresponding higher dimensional counterterms in the 4D theory 16 of 
zero-modes, ki = 0). Therefore, divergent terms q 4 R 2 /e can be avoided at one-loop provided 
that there is only one bulk propagator (see Table 1 for details when this can happen). 

With these considerations one concludes that the presence of q 4 R 2 /e and thus of the higher 
derivative counterterms is related to the multiplicity of the modes. Such operators are then due 
to the non-renormalisability of the models - initial supersymmetry cannot protect against their 
emergence as counterterms, regardless of the way supersymmetry is broken. Our analysis also 
shows that such operators are most relevant in models with low compactification scales, when 
their effects are less suppressed, see eas. tfTT?)) . (|41j) . 

Finally, note the presence of a u dependence of the coefficient of q R 2 /e in (|58|) . which is 
present in the case of brane- localised supersymmetry breaking (see Section 2.1). This signals 
some dependence between this supersymmetry breaking mechanism and the coefficient of the 
higher derivative operators. No such dependence appears for the discrete or continuous Scherk- 
Schwarz mechanism. 

4 Phenomenological implications: living with ghosts? 

The presence of higher derivative operators in the action of the scalar field 4>h has implications 

for phenomenology. Their investigation is however difficult since theories with higher derivative 

operators can bring in more fundamental problems, such as unitarity violation, non-locality 

effects, the presence of additional ghost fields, and for these reasons such theories were less 

15 see discussion in the last paragraph of this section. 
This is unlike the case of l/{R 2 e) divergences to which all modes including ki = contribute. 
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popular in the past (for some studies of such theories see for example |21|-|27p. Therefore the 
phenomenological considerations below should be taken with due care. 

In the presence of higher derivative operators, new fields (ghosts) are present. To see this, let 
us write the propagator of 4>h in the presence of the higher derivative operator found in eas.(|19|). 

<gH) : £ 4 = -£R 2 <p H n 2 (j>H + ■■■, where £ is an arbitrary constant (assumed positive for 
the convergence of the partition function); if <pH is a bulk field, we refer to its zero- mode only. 
The propagator of 4>h changes then into 



-£R 2 p 4 + p 2 



(1 -4£R 2 m 2 ) 



711 Z 



rrij 



with 



2£R 2 



l±(l-4£i? 2 m 2 



(62) 



(63) 



The second term in the rhs of (|62|) has the "wrong" sign, thus it signals the presence in the 
model of a ghost field of mass 777.4.. Here m 2 is the one- loop induced mass m;i (0) of 4>h plus the 
tree level contribution (if any), and mt\ is its value corrected by the higher derivative operator, 



but ignoring loop corrections 0(q 2 R 2 ) of eq.(J2J) at q 2 



■ m 



In function of the supersymmetry 



breaking mechanism which controls the coefficients C1.2 in (|16|) . one may have 111^(0) < and 
thus electroweak symmetry breaking (assuming no tree level mass is present). With m 2 < 0, 
£ > then rnz\_ > and mPz < 0, and the symmetry breaking may be maintained. 

In the following we require that > M 2 , i.e. the ghost mass m + is larger than the cutoff 
M* of our 5D theory, which for an effective theory approach as ours seems a natural requirement. 
We then study its implications for mz\. For this purpose we use a dimensional analysis |28l I29j 
to obtain perturbativity constraints on M*, eas. (|B-15|) . (|B-16|) . by requiring the effective gauge 
and Yukawa couplings be less than unity. The result is 

p-i 



12vr^ 1 



5 giR 1 



12vH 



2/(p-l) 



1 

R' 



(64) 



The first (second) bound is from gauge (Yukawa) interaction. Here p = 2 (p = 3) if the Higgs 
field is a brane (bulk) field; 5 is a factor equal to N for SU(N) gauge group, taking account of 
number of degrees of freedom present in loops j^Hi- Further, one imposes the first condition in 
the equation below, solves it for £ which is then used to evaluate the change in m 2 ^, to find: 

j3 sgn (m 2 ) 



Til 1 



7 V 



1 



7 2 2,4 



(65) 



sgn {m 2 ) 7 v 2 / (3 — 1 

Here M* ]U is upper bound of M* in v is the factor multiplying 1/R in M* )M , (3 = \m 2 \R 2 , 
sgn(x) is +1 (—1) for x > (x < 0), and A_ is the variation of m 2 _ relative to m 2 . One should 
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then take 7 > 1 but 7 ~ 0(1) may still satisfy > M 2 . We consider now the Yukawa 
correction only and take 17 5 = 3 for SU(3) C , m 2 < and f^t — 1 for which j3 ~ 0.1. This gives 
for 7 = 1, £ = 0.03(0.12) and A_ = -0.3% (-1.2%) for p = 2(3). Therefore the coefficient of 
the higher derivative counterterm must be very small and its correction to the scalar field mass 
is negligible. A value 7 = 1/4 would give £ = 0.13(0.53), A_ = -1.3% (-4.8%) for p = 2 (3), 
thus A_ is mildly changed. Note that this discussion ignored the additional terms iO{q 2 R 2 )) at 
q 2 = m 2 _ . For the implications on the physical mass of the scalar field and its dependence on the 
parameter £ one must actually evaluate the minimum of the one-loop potential computed in the 
presence of such operators, which is beyond the purpose of this work. Finally, if the condition 
m 2 ^ > M 2 is not satisfied, the presence of the ghost pole in the effective theory with the cutoff 
scale M* requires further theoretical investigation 18 . 

5 Conclusions. 

In this work we addressed the role that higher derivative operators play in the study of radiative 
corrections to the mass of the (Higgs) scalar field in 5D N=l supersymmetric models compactified 
on S\jZ<i- This is an important issue because it addresses how physics associated with compact 
dimensions decouples at low energies q 2 <C 1/R 2 , how these operators control the running of the 
scalar field mass across q 2 ~ 1 / R 2 , and the UV behaviour of the mass under the scaling of the 
momenta to q 2 3> 1/R 2 . Our technical results can be applied to a large class of orbifold models 
which investigate such one-loop effects from the bulk fields. 

The interactions considered were localised superpotentials and gauge interactions. Although 
the models are rather minimal, the multiplet structure and the interactions considered are 
generic in any realistic higher dimensional extensions of the SM or its supersymmetric versions. 
It was found that Yukawa interactions, unlike gauge corrections, can generate higher derivative 
counterterms to the scalar field mass even at the one-loop level. 

The work examined closely the relationship of the higher derivative operators with the way 
supersymmetry breaking was transmitted to the visible sector. The supersymmetry breaking 
scenarios addressed were local supersymmetry breaking by the F-term of a gauge singlet field 
localised at a hidden brane (fixed point) and also non-local breaking via Scherk-Schwarz (dis- 
crete/continuous) boundary conditions. In the case of local breaking, this effect is transmitted 
via radiative corrections to the scalar field, associated with the Higgs field. Such corrections, due 
to the compact dimension, are induced by the bulk fields which feel the supersymmetry breaking 
at the hidden brane. The second case, of non-local supersymmetry breaking, considered super- 
symmetry breaking by discrete and continuous Scherk-Schwarz twists of the bulk fields by using 

17 These correspond to the model with discrete Scherk-Schwarz supersymmetry breaking discussed above. 
18 for a discussion see [25) . 
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the SU(2)r symmetry of the multiplet content. 

Our analysis showed that in all these cases of supersymmetry breaking the emergence of 
higher derivative counterterms to the mass of the scalar field is very similar and has little or no 
dependence on the details of the breaking mechanism considered. This is the main result of the 
paper. At the technical level this means that the coefficient of these operators is independent 
of the boundary conditions for the bulk fields. As a result these operators seem to be a generic 
presence in orbifold compactifications and initial supersymmetry cannot protect against their 
presence, in some cases even at the one-loop level. This result is ultimately due to the non- 
renormalisable character of the theory and may also raise questions on the power of initial 5D 
supersymmetry in maintaining a mild UV behaviour in the compactified theory. 

The phenomenological implications of the presence of higher derivative operators in the 
action of the scalar field were briefly investigated. It was found that the requirement that the 
5D effective theory be weakly coupled together with the ghost field mass be larger than the 
5D effective cutoff lead to small corrections on the scalar field mass in the action. This can 
change dramatically when any of these constraints is relaxed. To evaluate the physical mass one 
is required to compute the one-loop corrected scalar potential in the presence in the action of 
higher derivative operators. Such calculation requires however a prior and more comprehensive 
study of the theories with higher derivative operators at the quantum level, which is beyond our 
purpose. 

Our study is also important because other approaches to compactification like string theory 
currently shed little light on such issues. The reason for this is that in string theory one computes 
the scalar potential (vacuum energy), derived at zero external momenta. Its second derivative 
giving a scalar mass cannot then recover a momentum dependence ("running") of the latter. 
We are thus confined to study these higher derivative operators in the framework of field theory 
orbifolds. The situation is very similar to the case of one-loop corrections to the gauge couplings, 
where again the role of higher derivative operators cannot be discussed in (on-shell) string loop 
corrections, but can be evaluated consistently in the context of field theory orbifolds |1UU111IT2| . 
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Appendix 

A Series of Integrals in the DR scheme. 

The functions used in eas. (jl6JI . I|17|). ifST))) are, up to a re-definition of e, 



Jp[c 1 ,c 2 ,c]= \ 



ni,rt2GZ ' 



dt 



-nt [c+ai(ni+ci) 2 +a2("2+C2) 2 ] 



, ai )2 > 0, c> 0; p= 1,2. (A-l) 



The expressions given below for J p , p = 1, 2 generalise results quoted in Appendix B of £I3 4 valid 
only for the case c/a\ 2 < 1 and which are included here for completeness. Using the method in 
Appendix B of |3U| one has that, if < c/a\ < 1 



Ji[ci,c 2 ,c] 



TIC 1 



In 47ra ie ^^i) + ^(- A ^) 



+ 27TM 



rl 
6 



+ \pKU > 



+ A* - (c/ai + 

rb+1/2] 



E ln |! 
meZ 



-27r 7(ni ) 



1P+1 



ai 



C[2p+l,l + A Cl ]+C[2p+l,l-A Cl ] , (A-2 



with u = vai/a2, 7 = 0.577216... If c/ai > 1, one has the same pole in e, but different finite part: 



Ji[ci,c 2 ,c] 



TTC 1 7TC 

H In 



l/rr ( "- J I - ^ In 1- 
mez 



-27T7(ni) 



+ 4 



c 

a 2 



— cos(2-7rniCi) K\ (^2-Khi{c/ a\) 2^ , 



rii>0 



(A-3) 



with the notation 



7(m) 



[z(ni)] 2 -ic 2 , 



z(ni) = c + ai(ni + c\) 2 . 



(A-4) 



Here £[2, a] is the Hurwitz Zeta function, (\z, a] = J2 n >o( n + a )~ z > Rez > 1, a ^ 0, —1, —2, • • •, 
and tp(x) = d/dx hiT[x]. Eqs. (|A-2|) . ljA-3|) depend on the fractional part of defined by 
A Ci = Ci — [cj] with 0< A Ci <1, [cj] G Z. Finally, is the modified Bessel function [3^ 
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dxx v - l e- bxP - ax - p = - 
p 



2p 



Kz(2Vab), Re(b), Re{a) > 



with 



(A-5) 
(A-6) 



which is strongly suppressed at large argument. 
One also finds that, if c <C 1 



Jl[ci,Cfc,C< 1] 

Ji[l/2,l/2,c< 1] 
J7i[l/2,0,c< 1] 
Ji[0,0,c< 1] 



7TC 1 



7TC 1 



:- - In 



7TC 1 



- - In 



1?1(C2 

— iuci|iit) _7,- nc 2 
(c2 — mci) ??(zm) 



^(l/2-m/2|m) ^ u/A 2 
ry(i u) 



ln(c/o2 + |C2 — UiCi| 



7TC 1 

\Ja\a 2 e 



In 



In 



$ 1 (-iu/2\iu) e _ W4 
(u/2) 7/(m) 

47T 2 |77(^)| 4 a 2 - x 



ln[u/2] : 



lnc, u = \ a\ja2 (A-7) 



used in the text. Above we used the Dedekind Eta function 



7y(r) 



3 7rrr/12 



2inT n 



), 



n>l 



v(~1/t) = ^i^v(r), ??(r + l) =e i7r / 12 7 ? (r), 
and the Jacobi Theta function $i 



(A-8) 



6i(z\t) = 2q 1 ^ sin(^) J] (1 - g n )(l - g n e 2i7rz )(l - g n e" 



g = e 



2J7TT 



n>l 



n i7rr(n+l/2) 2 (2n+l)iwz 



= -*£(-i) 

which has the properties 

<?i(i/|r + l) 
0i(i/ + 1|t) 

+ r | r) = — e 
0i(-i//t|-1/t) = e i7r/4 T 1/2 exp(^i/ 2 /r) ^i(^k) 



(A-9) 



e iw/4 $i{v\T) 



(A-10) 
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In the following we provide the results for J2 [c\ , C2 , c] whose calculation is almost identical. 
If < c/ai < 1 

ttV 1 vr 2 c 2 



J2[CI,C 2 ,C] 



2 ^/aia 2 e 2 ^/aio^ 



In 



+ ^3^-2Al(l + A^)-6- 

+ J] {(a2«(ni))3 Li 2 (e- 2 ^^)) + H Li 3 (e- 2 ^)) + C .cj 



4vrai e 7+V>(A cl )+v>(-A cl ) 
1 



- + A 2 

6 Cl 



c/ai + A, 



niGZ 

- w^rb+i/2] 



+ vr 3 / 2 a 2 U 3 £ 



p>i 



(p+2)! 



ai 



p+2 



(C[2p+1, 1+A ci ]+C[2p+1, 1-A C1 ]) .(A-ll) 



If c/ai > 1 the pole structure is similar: 



J 2 [ci,c 2 ,c] 



7T 2 C 2 1 7T 2 C 2 



:ln 



7T c e 



,7-3/2 



2 v /alai e 2 v /aia 2 
+ X! {(o2«(m))5 Li 2 (e- 2 ^)) + g Li 3 (e- 2 ^)) +c.c.| 



nieZ 



+ 4cu ^ cos(27rniCi) K 2 ^27rni(c/ai) 1 / 2 ^). 



fii>0 n i 



(A-12) 



For the simpler case c< 1, 



J 2 [ci,c 2 ,c< 1] 



7T 2 C 2 1 2 ol 
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--2A 2 (1-A C1 ) 2 
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niGZ 

with a ni = i c 2 + u |ni + ci|. Also 



2 ^/aia 2 e 

+ E {v^2"K+ci|Li 2 (e- 2 — i) + gLi 3 (e- 27rCT "i)+c. C .| (A-13) 



J 2 [0,0,c < 1] 



7T 2 C 2 1 7T 2 



+ — a 2 it 



2 ^/aia 2 e 45 
+ ^E {2vrn|n 1 |Li 2 ( e - 2 ^i)+Li 3 ( e - 2 ^i)} 



(A-14) 



ni6Z 
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with a ni = u \ n\\. Finally 

J 2 l/2,l/2,c«l = —a 2 u 3 

2 ^a\Q2 e ob(J 



niGZ 



+ V E { 27r u l n i + V2| Li 2 (-e- 2mT "i ) + Li 3 (-e- 2 ^i )} (A-15) 

with a ni = u\n\ + 1/2|. 



B Upper bounds on M* from Dimensional Analysis. 

We derive the upper bounds on the cutoff M*, used in the text, eas, ()64|) . For this purpose, we 
first outline the general dimensional analysis performed in f28 . (Section 3.2) and [22] (Section 3). 
We then apply it to our case to derive upper bounds on M* from perturbativity constraints. 
The action in D dimensions (for orbifolds with yi fixed points) has the form 

c = /W0) + E 6 °^y - w) fr) (b-i) 

i 

where <\> (ip) is a bulk (brane) field, respectively. One can assume canonical kinetic terms in (|B-1|I 
and then rescale these fields to their dimensionless counterparts <f>, if) and also the derivatives as 

tt*,y)= (l^-j <K x >y)> M*)=[iJfi) d = M*d, (B-2) 

where M* is the cutoff scale; Id is a suppression factor which accounts for angular integrations of 
loop corrections in D dimensions Id = (4:ir) D / 2 T(D /2) which grows rapidly with D, and 5 is the 
multiplicity of fields in loop diagrams for non-Abelian groups, for example 5 = N for SU(N), 
5 = 8 for 50(10). Using eq.JEH) in eq.CEQ) gives 

c = T ^c hulk ^) + J2^ D - 4 (y-y i )jj- 6 ^,A). (b-3) 

where £buik,£j only contain dimensionless couplings and fields. If all couplings in £buik, A are 
of order 1, all loops are of the same order of magnitude. The theory with £b u ik,£j remains 
weakly coupled if these dimensionless effective couplings remain smaller than unity [2H1 EH] . 
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Let us apply this result to (bulk) gauge interactions. The relation between the D dimensional 
gauge coupling and the 4D gauge coupling is 



Let us now consider a brane-localized interaction 



r£>-2\ \ ( -K/r2\\ 



(B-4) 



Vb-4 _ 1 

2 2 

ah at 

where Vd-a is the volume of extra dimensions. From the rescaled covariant derivative, 

V " ~ M. M. ~ M* WD \ l D /S ) ^' (B 5) 

we identify the dimensionless parameter corresponding to the gauge coupling as 

9 2 cS ^9l^j r <l. (B-6) 
The above condition for perturbativity imposes the bound on the cutoff as 



C F = J d 2 6W(<5>,y) + h.c. (B- 
where ^ (<&) is a brane [Zi even bulk) multiplet respectively. Under the rescaling 



*(x,y)=\-^-\ *(x,y), = ( j-J- 5 J c? 2 6> = M*d 2 9, 6 = "^|=> (B-9) 



M 4 

C F = —f- I d 2 6 W(4>, V) + h.c. (B-10) 



the brane action becomes 



Z 4 /<5 

In particular, for the superpotential for the Yukawa interaction 

W = X t H u QU c , with X t = V^ 4 U, t (B-ll) 

where p is the number of bulk fields present in the Yukawa interaction, and f±,t is the 4D 
coupling. Then, the redefined superpotential is given by 

W = A cff H u Q U c (B-12) 
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where 



v Z , , 3-p 

D-2 \ 2 / j\#2 \ ~~2~~ 



k/s M»-*y mi 

{ l D /5 ) \k/S 



< 1 



(B-13) 



This condition gives the bound on M* 



f2 yP 



p(D-4) 



(B-14) 



We now apply eqs. (|B-7j ). (|B-14|) to our 5D models, to derive bounds on the 5D cutoff M*. 
Ea.(|B-7|) gives 

M* < ^4- 4. (B-15) 



Eq. ()B-14|) becomes for the 5D case 



12vr 2 



16vr2 fl'^-V 



— 1 

p 1 



(B-16) 



where p = 2, 3. Eas. (|lTT5l) . (iB^Tnl) were used in the text, eas. (|64(l . 
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